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ABSTRACT 

We present a multi-wavelength study of embedded massive clusters in the nearby (3.9 Mpc) starburst 
galaxy NGC 4449 in an effort to uncover the earliest phases of massive cluster evolution. By combining 
high resolution imaging from the radio to the ultraviolet, we reveal these clusters in the process of 
emerging from their gaseous and dusty birth cocoons. We use Very Large Array (VLA) observations at 
centimeter wavelengths to identify young clusters surrounded by ultra-dense Hll regions, detectable via 
their production of thermal free-free radio continuum. Ultraviolet, optical and infrared observations 
are obtained from the Hubble and Spitzer Space Telescope archives for comparison. We detect 39 
compact radio sources towards NGC 4449 at 3.6 cm using the highest resolution (1'.'3) and sensitivity 
(RMS ^ 12 /iJy) VLA image of the galaxy to date. We reliably identify 13 thermal radio sources 
and their physical properties are derived using both nebular emission from the Hll regions and SED 
fitting to the stellar continuum. These radio detected clusters have ages < 5 Myr and stellar masses 
of order 10* Mq. The measured extinctions are quite low: 12 of the 13 thermal radio sources have 
Ay < 1.5, while the most obscured source has Ay w 4.3. By combining results from the nebular and 
stellar emission, we find an /-band excess that is anti-correlated with cluster age and an apparent 
mass-age correlation. Additionally, we find evidence that local processes such as supernovae and 
stellar winds are likely playing an important role in triggering the current bursts of star formation 
within NGC 4449. 

Subject headings: galaxies: dwarf - galaxies: individual (NGC 4449) - galaxies: irregular - galaxies: 
starburst - galaxies: star clusters 



1. INTRODUCTION 

Young massive star clusters represent a fundamental 
mode of star formation throughout the universe, both lo- 
cally and at high rcdshift. The most massive and dense 
clusters are consistent with being young analogues of an- 
cient globular clusters like those in our own Milky Way 
(e.g. Ho & Filippenko 1996a,b). The prevalence of glob- 
ular clusters around massive galaxies today suggests that 
these extreme star clusters were forming ubiquitously in 
the early universe as galaxies began to merge and coa- 
lesce hierarchically. 

Young massive clusters (YMCs) continue to play an 
important role in the evolution of the universe. YMCs 
and their constituent massive stars ionize the interstel- 
lar medium, power the infrared radiation of dust, pro- 
duce heavy elements, and trigger further star formation 
via supernova explosions. In addition, the most mas- 
sive YMCs, known as "super star clusters" (see Whit- 
more 2003; O'Connell 2004, and reference therein for 
overviews), are packed with hundreds to thousands of 
O and B stars. These extreme clusters can drastically 
alter the morphology of their host galaxy when the mas- 
sive stars collectively explode at the end of their lives, 
expelling huge amounts of gas and blowing galactic-scale 
superbubbles (e.g. Tenorio-Tagle et al. 2007; Marlowe et 
al. 1995). 

Despite the importance of massive stars and massive 
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clusters, their formation and earliest evolutionary stages 
are not well understood (Tan 2005; Zinnecker & Yorke 
2007). We do know, however, that nearly all massive 
stars form in clusters and that YMCs host large numbers 
of massive stars (Hunter 1999; Clark et al. 2005; de Wit 
et al. 2005). Therefore, in addition to being interesting in 
their own right, YMCs provide important laboratories for 
studying the clustered mode of massive star formation, a 
crucial part of understanding massive star formation in 
general. 

Interest in massive extragalactic star clusters was 
sparked by the discovery of a population of YMCs in the 
peculiar galaxy NGC 1275 by Holtzmann et al. (1991) 
using the Hubble Space Telescope (HST). Since then, 
young massive clusters have been found in a large num- 
ber of galaxy systems, primarily at optical wavelengths 
with HST (for reviews, see Whitmore 2003; Larsen 2006). 
More recently, ultra-young (< few Myr) massive clusters, 
still deeply embedded in their birth material, have been 
discovered using infrared and radio imaging (e.g. Kobul- 
nicky & Johnson 1999; Turner et al. 2000; Neff & Ulves- 
tad 2000; Johnson et al. 2001; Beck et al. 2002; Johnson 
et al. 2003; Johnson & Kobulnicky 2003; Johnson et al. 
2004; Tsai et al. 2006). By combining observations span- 
ning multiple wavelength regimes, we aim to improve the 
current understanding of massive cluster (and massive 
star) formation and early evolution. 

Here we present our study of embedded massive 
clusters in the nearby starburst galaxy NGC 4449. 
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TABLE 1 
VLA Observations of NGC 4449 



Wavelength Weighting Orig. Synthesized Beani^ 

(cm) (robust value) (arcscc X arcsec) 



P.A. RMS Noise 

(degrees) (fJ-Jy bcani^^) 



6.0 


5 


0.70 X 0.65 


83 


36 


3.6 


1 


1.27 X 1.14 


-73 


12 


1.3 


5 


1.29 X 1.04 


-86 


25 



^The data have been convolved to identical beams of 1'.'3 X 1'.'3 to better facilitate the 
comparison between wravelengths. 



NGC 4449 is a barred Magellanic-type irregular galaxy 
(de Vaucouleurs et al. 1991) with star formation occur- 
ring throughout the galaxy at a rate almost twice that 
of the Large Magellanic Cloud (Thronson et al. 1987; 
Hunter et al. 1999). NGC 4449 has almost certainly 
undergone an external perturbation, as is evident in Hi 
studies revealing extended streamers wrapping around 
the galaxy and counter-rotating (inner and outer) gas 
systems (Hunter 1998; Theis & Kohle 2001, and refer- 
ences therein). The proximity of NGC 4449 is advan- 
tageous for extragalactic star cluster research since rel- 
atively small spatial scales can be probed. The work of 
Gelatt et al. (2001) has already revealed ^ 60 optically- 
selected compact star clusters. We adopt a distance of 
3.9 Mpc, consistent with recent work by Annibali et al. 
(2007), and thus 1"=19 pc. 

Using high-resolution radio imaging, we identify the 
youngest embedded star clusters in NGC 4449. The ob- 
served ultra-dense Hll regions are produced when hot 
massive stars within the clusters ionize the surrounding 
gas and are identified as compact free-free (thermal) ra- 
dio sources. A previous catalog of radio detected Hll 
regions in NGC 4449 was presented by Israel (1980) us- 
ing an NRAO 3-element image at 11 cm (2.7 GHz) from 
Seaquist & Bignell (1978) with a beam size of 10'.'8 x 6'.'5 
and an RMS noise of 0.4-0.6 mJy beam~^. Israel (1980) 
derived a list of eight Hll regions that had counterparts 
in the optical Ha data of Crillon & Monnet (1969). The 
optical sizes were 10"-20" and the 11 cm flux densities 
were 1-10 mJy. Israel & Kennicutt (1980) also calculated 
the radio to Ha derived extinction for seven objects and 
compared these values with the Balmer decrement de- 
rived extinctions for a single object (CM 39). 

Compared to the NRAO 3-element image at 11 cm 
used by Israel (1980), the current Very Large Array 
(VLA) survey at 3.6 cm has a sensitivity that is im- 
proved by a factor of ^ 40 and a beam solid angle that 
is smaller by a factor of ^ 40. In the present work, 
we combine high-quality VLA data and archival images 
from the Hubble and Spitzer Space Telescopes with stel- 
lar population synthesis models in order to discern the 
physical properties of the radio-detected star clusters in 
NGC 4449 and gain insight into the earliest evolutionary 
phases experienced by massive star clusters. 

2. OBSERVATIONS 

2.1. Radio Imaging with the VLA 

Radio observations of NGC 4449 were obtained with 
the Very Large Array^ from April 2001 to November 

^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 



2002. Observations at C-band (6 cm, 5 GHz) were 
carried out with the A-array, observations at X-band 
(3.6 cm, 8 GHz) were carried out with both the A- and 
C-array, and observations at K-band (1.3 cm, 23 GHz) 
were carried out with the C-array. The sources 3C48 and 
3C286 were used as flux density calibrators, and we esti- 
mate that the absolute flux density calibration is < 5%, 
based on the scatter in the VLA Flux Density Calibrator 
database. The high-frequency observations at K-band 
utilized fast-switching to a nearby phase calibrator, dis- 
paced by ~ 5°, in order to mitigate the effect of the 
atmosphere. 

These data were reduced and calibrated using the As- 
tronomical Image Processing System (AIPS). In order to 
better match the largest spatial scale to which the ob- 
servations at each wavelength is sensitive, antenna spac- 
ings of < 20 kA were not used. Additionally, in order 
to achieve synthesized beams of approximately the same 
size at each wavelength, the three data sets were im- 
aged with different weighting of the uv-coverage, ranging 
from purely natural weighting at 6 cm to slightly uniform 
weighting at 1.3 cm. The data were convolved to identi- 
cal beams of 1'.'3 x 1'.'3 in the imaging process to better 
facilitate comparison between wavelengths. Despite our 
best efforts, it should be noted that it is not possible 
to precisely match the uv-coverage at each wavelength, 
and there are slight variations in sensitivity to different 
spatial scales. After imaging the data, a primary beam 
correction was applied at each wavelength. The resulting 
parameters are listed in Table 1. 

We identify 39 compact radio sources towards 
NGC 4449 which are shown in Figure 1 and labeled in 
order of increasing right ascension. Source identification 
is based on a minimum 3a detection (locally) in the 3.6 
cm image, which has the highest sensitivity of the three 
radio images. Figure 2 shows an archival HST F550M 
(^ t^-band) image in yellow overlaid with our VLA 3.6 
cm image in blue. 

2.2. Ultraviolet, Optical and Near-infrared Imaging with 

HST 

A wealth of data on NGC 4449 is contained in the 
HSI^ archive and we have obtained broad and narrow- 
band images of the galaxy spanning ultraviolet to near- 
infrared wavelengths. These observations were taken 

merit by Associated Universities, Inc. 

^ Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained [from the Data Archive] at the Space 
Telescope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under NASA con- 
tract NAS 5-26555. These observations are associated with pro- 
gram #6716, 7919, 10522 & 10585. 
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Fig. 1. — VLA 3.6 cm image of NGC 4449. The radio sources have been labeled in order of increasing right ascension. This image has a 
resolution of I'.'S, a sensitivity of ^ 12 /ijy (RMS), and has been corrected for the primary beam. 



over a period of several years with the Advanced Cam- 
era for Surveys (ACS), the Wide Field and Planetary 
Camera 2 (WFPC2), and the Near Infrared Camera and 
Multiobject Spectrometer (NICMOS). The /fS'T obser- 
vations have spatial resolutions ranging from ^ 0'.'l-0'.'2. 
A summary of the if5'T observations used in this project 
is given in Table 2. 

The ACS Wide Field Channel (WFC) has a field of 
view of 202" X 202" and a plate scale of - 0'.'05 pixer^. 
The ACS pipeline employs the MultiDrizzle (Koekemoer 



et al. 2002) software to produce cosmic ray cleaned, 
combined images and we have obtained these data on 
NGC 4449 through the F435W, F550M, F814W, F658N, 
and F660N filters. These observations provide coverage 
of all of the locations of the radio sources in NGC 4449 
identified in Figure 1. 

WFPC2 is composed of four CCDs: three Wide Field 
(WF) chips and the Planetary Camera (PC). The plate 
scales of the PC and WF chips are ^ 0"05 pixel^^ and 
^ O'.'l pixel^^, respectively. The WF chips provide an 
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Fig. 2. — An HST/ ACS F550M image (yellow) overlaid with our VLA 3.6 cm image (blue) highlighting sites of current star formation 
in NGC 4449. The field of view of this image is ~ 2f6 X 3.'0. 



L-shaped field of view of 150" x 150" and the PC fills in 
another 35" x 35". All but two of the radio sources (6 
& 29) in NGC 4449 have WFPC2 coverage through the 
F170W and F336W filters. Two exposures were taken 
through each of these filters and we combined the cali- 
brated images, chip by chip, and rejected cosmic rays. 

A snapshot NICMOS observation through the F160W 
filter was obtained with the NIC 3 Camera and we have 
obtained this calibrated image from the archive. This 
image, with a field of view of 51"x 51" and plate scale 



of ~ 0"2 pixel^^, only covers 7 radio sources (15, 16, 17, 
19, 23, 24 & 25) in the nuclear region of the galaxy. 

We register all of the HST images to match the as- 
trometry of the VLA images, which is accurate to ~ 0"1. 
Coincident sources (typically 5-10) are found in pairs of 
images and the mean pixel offsets are found between the 
two images. We begin by registering the F658N (Ha) 
image to the 3.6 cm radio continuum image since both 
observations trace Hll regions. The other HST images 
can then be registered with respect to this or subsequent 
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TABLE 2 
Archival HST Observations of NGC 4449 



Filter 


Instrument 


Deseription 


Date Observed 


Proposal ID 


PI 




Broad and Medium-band Filters 


F170W 
F336W 
F435W 
F550M 
F814W 
F160W 


WFPC2 

WFPC2 

ACS/WFC 

ACS/WFC 

ACS/WFC 

NIC 3 


UV 

WFPC2 U 

Johnson B 

Narrow V 

Broad I 

H 


1997 Jul 28, 1998 Jan 9 

1997 Jul 28, 1998 Jan 9 

2005 Nov 10-11 

2005 Nov 18 

2005 Nov 10-11 

1998 June 13 


6716 
6716 
10585 
10522 
10585 
7919 


T. Steeher 
T. Steeher 

A. Aloisi 
D. Calzetti 

A. Aloisi 
W. Sparks 


Narrovir-band Filters 


F658N 
F660N 


ACS/WFC 
ACS/WFC 


Ha -1- [NII]A6584 
[NII]A6584 


2005 Nov 17 
2005 Nov 18 


10585 
10522 


A. Aloisi 
D. Calzetti 



registered images. The HST images are aligned with re- 
spect to the VLA images to a precision of ^ 0'.'2. 

2.3. Infrared Imaging with Spitzer 

We have obtained infrared images (3.6-24 /im) of 
NGC 4449 from the Spitzer Space Telescope archive. We 
primarily use these data in this paper to qualitatively in- 
vestigate the infrared properties of NGC 4449 since these 
images suffer from significantly lower resolution (~ 3-6") 
than the VLA and HST images. However, despite their 
limitations, the infrared data are informative, providing 
the general IR morphology of NGC 4449 and highlighting 
regions associated with warm dust. 



2.4. 



A Multi-wavelength View of NGC 4449 



The high quality data used in this study of embedded 
massive clusters allows for a detailed inspection of the 
starburst morphology of NGC 4449. Figure 3 shows a 
selection of images of the galaxy at multiple wavelengths. 
"Close-up" 3-color images (3.6 cm, Ha and 1^-band) are 
shown in Figure 4 with the sources labeled according to 
Figure 1. 

The similarities between the 3.6 cm and Ha observa- 
tions are immediately apparent in Figure 3. This like- 
ness is not surprising since both observations trace dense 
ionized gas. The 3.6 cm and 24 /xm images are also strik- 
ingly similar, indicating that the presence of warm dust 
is almost exclusively associated with the ultra-dense Hll 
regions in the galaxy. 

The circles in Figure 3 are centered on two example 
radio sources: one with and one without counterparts in 
the other wavebands. Source 30 (to the NE) has both 
optical and infrared counterparts, suggesting this stellar 
cluster is in the process of emerging from its birth mate- 
rial, all the while ionizing the surrounding gas and heat- 
ing the dust. Source 28 on the other hand, does not have 
detectable counterparts in any of the other wavebands. 
We will return to the radio sources without optical or 
infrared counterparts in §5.4. 

3. MULTI- WAVELENGTH FLUX MEASUREMENTS OF THE 
RADIO SOURCES 

3.1. Methodology 

In order to study the emergence of young massive star 
clusters, we first use radio continuum imaging to identify 



potential clusters surrounded by ultra-dense Hll regions 
(see Fig. 5 for some examples). We then measure the 
ultraviolet to radio fiux densities of the radio-selected 
sources (in the VLA and HST images) using a new pho- 
tometry code that we developed (SURPHOT), allowing 
for consistent irregular apertures and background annuli 
across multiple wavebands. SURPHOT is first used to 
detect a specified contour level in a selected region of the 
3.6 cm image which is used as a reference. A source is 
selected and the contour may be used for the aperture in 
all of the images. Alternatively, an ellipse, circle or free- 
form aperture can be used if appropriate (e.g. nearby 
yet distinct sources that are contained within the same 
contour). A background annulus is constructed by radi- 
ally expanding the chosen aperture by specified factors 
(typically 1.75x and 2.5x). The pixel coordinates of 
the aperture and background annulus are transformed 
between images using their FITS World Coordinate Sys- 
tems (WCS) which relate image pixels to sky coordinates. 
In this way, flux densities are computed in all (11) of 
the VLA and HST images using identical apertures and 
background annuli. (A more detailed description of SUR- 
PHOT is provided in the Appendix.) The resulting flux 
densities are used in conjunction with stellar population 
synthesis models to derive physical properties such as 
extinctions, ages, and masses. 

Our choice to use radio-selected apertures is both sci- 
entiflcally and practically motivated. The focus of this 
work is on extremely young massive clusters and cen- 
timeter radio emission is a good indicator of youth. In 
addition, the resolution of the VLA images (1'.'3) is larger 
than the resolutions of the iJ6'T images (~ 0'.'l-0'.'2). We 
note that using a larger aperture than the intrinsic size of 
a stellar cluster (in the HST data,) allows for the possibil- 
ity of contamination from an underlying or surrounding 
stellar population within the apertures. However, this 
effect is not likely to be very signiflcant since the stel- 
lar light is domiated by the clusters of interest and we 
perform background subraction in the photometry. 

3.2. The Radio Sources 

Thermal radio sources are recognized as having flat 
or positive radio spectral energy distributions {Si, ~ 
i^^jCt > 0), indicative of free- free emission from dense 
ionized gas associated with young massive star clusters. 
Non-thermal sources have decreasing radio SEDs charac- 
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Fig. 3. — Multi-wavelength images of NGC 4449 (FOV ^ 2.'5 X 2.'5). The following observations are displayed: 3.6 cm (radio continuum), 
F658N (Ha), F550M (F-band), 4.5 /im (stellar continuum & warm dust), 8 /im (warm dust & polycyclic aromatic hydrocarbons (PAHs)) 
and 24 ^m (warm dust). The circles (r = 3") surround two radio sources: Source 30, to the NE, has counterparts in all of the other 
wavebands and has therefore already begun to emerge from its gaseous and dusty birth cocoon. Source 28, on the other hand, does not 
have counterparts in any of the other images. The nature of this source is unclear. 



teristic of supernova remnants and (background) active 
galactic nuclei. Here we are concerned with the ther- 
mal sources in NGC 4449, which must be disentangled 
from the non-thermal sources by calculating the spec- 
tral indices, 
densities. 



^3. Gem 



and af'3^™, using the measured flux 



We use the 3.6 cm radio image as our reference for 
both source identification and as the reference image in 
SURPHOT. The 1.3 cm image is somewhat less likely to 
suffer from non-thermal contamination or self-absorption 
but the sensitivity is lower than the 3.6 cm image by a 
factor of ^ 2 at best (the image quality is degraded far- 
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Fig. 4. — Close-up 3-color images of radio-detected sources in 
NGC 4449: 3.6 cm (green), Hq (red), K-band (blue). At the dis- 
tance of this galaxy (3.9 Mpc), 1"= 19 pc. 



Fig. 5. — Examples of radio-detected star clusters in NGC 4449. 
HST/ ACS 3-color images ('^B, ~ V, ~/) are shown with 3cr radio 
(3.6 cm) contours overlaid. 

ther from the center of the 1.3 cm image). Source selec- 
tion is based on a 3cr detection threshold, where a is the 
local RMS noise (see Fig. 1). We estimate the minimum 
uncertainty in the measured radio flux densities as this 
local RMS noise. In many cases, however, the dominant 
source of uncertainty comes from the choice of aperture. 
For faint sources, estimates of the background level can 
also lead to large uncertainties. We therefore make mul- 
tiple flux measurements for each source, varying these 
parameters. The fluctuations in these measurements are 
accounted for in our final reported values of flux densi- 
ties and spectral indices, which are given in Table 3 along 
with the coordinates of the 39 radio sources identified in 
NGC 4449. In the cases of weak or non-detections in 
the 6.0 and 1.3 cm images, 3(T detection limits are given 
for the fiux densities. Corresponding limits are given for 
the spectral indices. Sources are designated as either 
(likely) thermal, non-thermal, or mixed/uncertain based 
on their spectral indices as described above. Of the 39 
radio sources, 13 are thermal, 22 are uncertain or com- 
posed of mixed sources and 4 are non-thermal. 

Of the 4 definitive non-thermal sources, 3 appear to be 
associated with the galaxy (Sources 13, 14, & 39) while 
Source 6 is most likely a background AGN with no optical 
counterpart. Source 14 is a well-studied, young (60-200 
years) luminous SNR (J1228-f 441) (see Lacey et al. 2007, 
and references therein). Lacey et al. (2007) have given 
a modern summary of the 30 year history of the radio 
hght curve of J1228-I-441. Source 13 may also be a SNR, 
although it appears to be at least partially embedded in 
an Hll region. Source 39 is quite intriguing. The radio 
emission is partially associated with an obvious point- like 
optical counterpart (probably a SNR) within a larger bi- 
conical shaped Hll region (see Fig. 4). The non-thermal 
sources within NGC 4449 are interesting in their own 
right, but they are excluded from further analysis in this 
work. 

3.3. Photometry of the Radio Sources in the HST Data 

We use the 3.6 cm radio map as the reference image 
in SURPHOT (see Appendix) to obtain photometry of 
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TABLE 3 
Radio Sources in NGC 4449: Coordinates. Flux Densities and Spectral Indices 



Source 


R.A. 


Decl. 


Se.Ocm 


Ss.&cin 


Si. 3cm 


^6. 0cm 


^3. 6cm 
*^1.3cm 


Designation 




(J2000.0) 


(J2000.0) 


(MJy) 


(/^Jy) 


(MJy) 








1 


12 28 08.86 


+44 05 09.0 


< 90 


80(20) 


< 210 


> -0.3 


< 1.1 


uncertain 


2 


12 28 09.12 


+44 05 20.9 


< 110 


120(20) 


< 150 


> 0.2 


< 0.3 


likely thermal 


3 


12 28 09.37 


+44 05 20.2 


330(40) 


280(30) 


280(50) 


-0.3(0.3) 


0.0(0.2) 


likely thermal 


4 


12 28 09.44 


+44 05 16.3 


230(40) 


270(30) 


330(60) 


0.3(0.4) 


0.2(0.2) 


thermal 


5 


12 28 09.70 


+44 05 19.8 


< 110 


130(30) 


< 140 


> 0.3 


< 0.1 


likely thermal 


6 


12 28 09.70 


+44 05 43.5 


410(30) 


250(10) 


140(30) 


-0.9(0.2) 


-0.6(0.2) 


non-thermal 


7 


12 28 10.01 


+44 05 24.7 


< 120 


130(20) 


170(40) 


> 0.2 


0.3(0.3) 


thermal 


8 


12 28 10.38 


+44 06 46.1 


< 140 


150(20) 


< 120 


> 0.2 


< -0.2 


uncertain 


9 


12 28 10.56 


+44 05 31.5 


< 120 


70(20) 


< 110 


> -0.8 


< 0.4 


uncertain 


10 


12 28 10.66 


+44 06 49.3 


< 140 


160(20) 


< 130 


> 0.3 


< -0.2 


uncertain 


11 


12 28 10.81 


+44 05 34.1 


160(40) 


230(30) 


220(40) 


0.6(0.5) 


0.0(0.2) 


thermal 


12 


12 28 10.82 


+44 05 36.8 


210(50) 


130(20) 


150(40) 


-0.8(0.5) 


0.1(0.3) 


mixed 


13 


12 28 10.87 


+44 06 51.5 


360(50) 


110(20) 


< 140 


-2.1(0.4) 


< 0.2 


non-thermal 


14 


12 28 10.94 


+44 06 48.5 


5,920(100) 


2,120(30) 


900(50) 


-1.8(0.0) 


-0.9(0.1) 


non-thermal'' 


15 


12 28 11.30 


+44 05 38.1 


2,270(60) 


2,250(60) 


1,510(40) 


0.0(0.1) 


-0.4(0.0) 


mixed 


16 


12 28 11.85 


+44 05 41.0 


170(50) 


110(20) 


100(30) 


-0.7(0.6) 


-0.1(0.4) 


mixed 


17 


12 28 12.40 


+44 05 47.1 


< 120 


110(20) 


< 80 


> -0.1 


< -0.3 


uncertain 


18 


12 28 12.63 


+44 05 03.7 


< 90 


150(10) 


< 160 


> 0.8 


< 0.1 


likely thermal 


19 


12 28 12.79 


+44 06 12.0 


< 110 


60(10) 


< 70 


>-1.0 


< 0.1 


uncertain 


20 


12 28 12.82 


+44 05 04.4 


< 90 


40(10) 


< 150 


> -1.4 


< 1.3 


uncertain 


21 


12 28 12.91 


+44 04 59.0 


< 100 


80(10) 


< 180 


> -0.3 


< 0.8 


uncertain 


22 


12 28 12.99 


+44 06 56.3 


150(40) 


160(20) 


180(50) 


0.1(0.5) 


0.1(0.3) 


thermal 


23 


12 28 13.08 


+44 05 42.7 


500(60) 


470(20) 


480(50) 


-0.1(0.2) 


0.0(0.1) 


thermal 


24 


12 28 13.74 


+44 06 09.6 


< 110 


100(10) 


< 60 


> 0.0 


< -0.5 


uncertain 


25 


12 28 13.76 


+44 06 12.3 


< 110 


90(10) 


< 70 


> -0.3 


< -0.3 


uncertain 


26 


12 28 13.86 


+44 07 10.4 


1,970(40) 


1,800(40) 


1,360(90) 


-0.2(0.1) 


-0.3(0.1) 


mixed 


27 


12 28 14.83 


+44 07 10.0 


290(40) 


300(20) 


340(90) 


0.1(0.3) 


0.1(0.3) 


thermal 


28 


12 28 14.88 


+44 06 23.2 


< 110 


80(10) 


< 70 


> -0.5 


< -0.1 


uncertain 


29 


12 28 15.09 


+44 07 10.5 


< 110 


160(30) 


< 240 


> 0.6 


< 0.4 


likely thermal 


30 


12 28 16.02 


+44 06 29.3 


130(80) 


190(20) 


160(50) 


0.7(1.1) 


-0.2(0.3) 


likely thermal 


31 


12 28 16.13 


+44 06 42.6 


< 100 


50(10) 


< 120 


>-1.2 


< 0.9 


uncertain 


32 


12 28 16.14 


+44 06 44.5 


< 100 


60(20) 


< 120 


> -0.9 


< 0.7 


uncertain 


33 


12 28 16.26 


+44 06 40.8 


< 100 


120(20) 


< 130 


> 0.3 


< 0.0 


likely thermal 


34 


12 28 16.32 


+44 06 38.3 


< 100 


50(10) 


< 120 


> -1.2 


< 0.8 


uncertain 


35 


12 28 16.41 


+44 06 43.6 


< 100 


50(10) 


< 130 


>-1.3 


< 0.9 


uncertain 


36 


12 28 16.49 


+44 06 44.7 


< 100 


40(10) 


< 130 


>-1.5 


< 1.1 


uncertain 


37 


12 28 17.69 


+44 06 29.1 


< 100 


240(20) 


< 150 


> 1.6 


< -0.5 


uncertain 


38 


12 28 18.13 


+44 06 35.6 


< 100 


90(20) 


< 180 


> -0.3 


< 0.7 


uncertain 


39 


12 28 20.14 


+44 06 16.2 


220(30) 


140(20) 


< 320 


-0.8(0.4) 


< 0.9 


non-thermal 


Note. — Units of right ascension are hours, minutes, and seconds. 


and units of declination are degrees, arcminutes, 


and arcseconds. The uncertainties in the right 


ascension and declination are 0.01 seconds and 0.1 arcseconds, respectively. 



Source designations are based on the radio spectral indices as described in the text. 
''Source 14 is a well-studied supernova remnant in NGC 4449 (e.g. Lacey et al. 2007). 



the radio sources in the HST data. HST system mag- 
nitudes with Vega zeropoints are reported in Table 4. 
Error estimates include contributions from Poisson noise 
and uncertainties from varying the apertures (see §3.2). 

4. PROPERTIES OF THE EMERGING STAR CLUSTERS IN 
NGC 4449 

4.1. The Models 

In order to estimate the physical properties of the ra- 
dio selected sources in NGC 4449, we compare our data 
to the latest starburst99 population synthesis models 
(Version 5.1) of Leitherer et al. (1999). We have run 
simulations at two metallicities, Z=0.004 and Z=0.008, 
as appropriate for NGC 4449 (Lequeux et al. 1979). We 
adopt an instantaneous burst of 10^ Mq with a Kroupa 
IMF, the Geneva evolutionary tracks with high mass loss, 
and the Pauldrach/Hillier atmospheres. Ages, masses, 
and extinctions are estimated using two distinct meth- 
ods: 1) comparing the nebular emission from the Hll 
regions to the model output (§4.2); and 2) fitting model 
spectral energy distributions (SEDs) to the broad-band 



FS-T flux densities (§4.3). 

4.2. Results from the Nebular Emission 
4.2.1. Ages 

We estimate ages of the radio sources in NGC 4449 
by comparing their measured Ha equivalent widths with 
the model predictions as a function of age. For an instan- 
taneous burst of star formation, the equivalent width of 
Ha emission is a reliable age indicator for young clus- 
ters since it measures the ratio of the ionizing flux (from 
massive stars) to the total continuum flux (e.g., Alonso- 
Herrero et al. 1996). This ratio is strongly dependent on 
age for stellar populations < 10 Myr, so long as the most 
massive stars have begun to die, thereby decreasing the 
ionizing flux. In addition, if the extinction is the same 
to the stars and the gas, Ha equivalent width will be 
independent of extinction. 

We calculate Ha equivalent widths, VFhq, using Equa- 
tions 1 & 2 below: 



Source 
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TABLE 4 
HST Vega Magnitudes of the Radio Sources in NGC 4449 



F170W 



F336W 



F435W 



F550M 



F814W 



F160W 



F658N 



F660N 



1 


>18.40 


19.93(0.35) 


21.47(0.08) 


21.37(0.09) 


20.05(0.05) 






16.72(0.02) 


18.26(0.07) 


2 


>18.00 


>20.53 


>22.66 


>22.51 


>21.71 






18.96(0.19) 


>21.29 


3 


15.86(0.08) 


16.57(0.07) 


17.65(0.17) 


17.76(0.12) 


17.18(0.15) 






14.30(0.03) 


15.76(0.04) 


4 


14.59(0.04) 


15.28(0.03) 


16.74(0.02) 


16.67(0.03) 


16.64(0.03) 






14.15(0.04) 


15.38(0.04) 


5 


16.82(0.13) 


17.50(0.05) 


18.71(0.05) 


18.88(0.05) 


18.43(0.08) 






16.01(0.09) 


16.96(0.08) 


6 






21.02(0.31) 


19.88(0.10) 


19.92(0.42) 






20.26(0.40) 


19.95(0.10) 


7 


16.59(0.27) 


17.92(0.12) 


19.42(0.08) 


19.50(0.03) 


19.04(0.14) 






15.82(0.05) 


17.40(0.05) 


8 


16.94(0.11) 


17.39(0.05) 


18.95(0.04) 


19.12(0.07) 


18.60(0.04) 






15.14(0.07) 


16.93(0.08) 


9 


>18.16 


17.95(0.03) 


19.34(0.08) 


19.05(0.08) 


18.30(0.06) 






17.38(0.07) 


18.06(0.08) 


10 


17.71(0.20) 


17.74(0.09) 


19.28(0.27) 


18.83(0.04) 


18.66(0.23) 






15.78(0.08) 


17.35(0.07) 


11 


>17.79 


17.09(0.26) 


17.73(0.09) 


18.04(0.26) 


16.92(0.10) 






15.37(0.08) 


16.55(0.17) 


12 


>17.73 


18.65(3.05) 


18.23(0.08) 


17.86(0.27) 


17.23(0.06) 






15.70(0.16) 


16.78(0.26) 


13 


>18.73 


18.90(0.07) 


20.38(0.11) 


19.95(0.19) 


19.70(0.14) 






16.04(0.08) 


17.63(0.07) 


14 


16.40(0.15) 


17.18(0.12) 


18.31(0.03) 


18.67(0.15) 


17.81(0.02) 






16.61(0.20) 


17.51(0.11) 


15 


13.48(0.08) 


13.81(0.03) 


14.64(0.09) 


14.29(0.05) 


13.62(0.07) 


12.06(0.04) 


12.70(0.05) 


13.37(0.04) 


16 


>18.09 


18.35(0.10) 


19.40(0.16) 


19.29(0.05) 


18.26(0.11) 


16.89(0.14) 


16.88(0.03) 


17.91(0.17) 


17 


17.54(0.24) 


17.72(0.07) 


18.48(0.15) 


18.31(0.08) 


17.63(0.14) 


16.23(0.10) 


15.40(0.04) 


16.85(0.06) 


18 


16.11(0.04) 


17.01(0.02) 


18.65(0.03) 


18.91(0.01) 


18.41(0.05) 




14.91(0.02) 


16.77(0.02) 


19 


>18.58 


19.34(0.13) 


21.00(0.39) 


20.43(0.25) 


20.70(0.49) 


19.01(0.36) 


17.23(0.18) 


18.52(0.15) 


20 


>18.60 


20.56(0.34) 


>23.42 


>23.23 


>22.62 




18.60(0.08) 


19.20(0.27) 


21 


15.94(0.04) 


17.03(0.02) 


18.56(0.02) 


18.89(0.02) 


18.47(0.05) 




15.22(0.04) 


16.94(0.03) 


22 


15.59(0.03) 


16.49(0.01) 


18.08(0.02) 


17.88(0.01) 


17.91(0.01) 




15.19(0.02) 


16.66(0.01) 


23 


16.46(0.15) 


16.54(0.04) 


18.22(0.06) 


18.68(0.03) 


17.57(0.04) 


16.00(0.09) 


14.46(0.05) 


15.70(0.04) 


24 


17.90(0.23) 


18.71(0.15) 


19.72(0.13) 


19.38(0.17) 


18.92(0.21) 


17.05(0.21) 


17.94(0.14) 


17.95(0.10) 


25 


16.71(0.10) 


17.48(0.06) 


19.25(0.09) 


19.00(0.10) 


18.75(0.10) 


17.49(0.11) 


15.75(0.06) 


17.02(0.05) 


26 


14.03(0.01) 


14.78(0.03) 


16.20(0.00) 


16.44(0.02) 


15.84(0.01) 






12.65(0.02) 


14.41(0.02) 


27 


16.44(0.06) 


17.04(0.03) 


18.29(0.08) 


18.44(0.02) 


17.87(0.09) 






14.70(0.03) 


16.32(0.04) 


28 


>18.84 


>21.34 


22.03(0.07) 


21.18(0.09) 


20.51(0.16) 






20.92(0.17) 


>21.67 


29 






19.59(0.43) 


19.57(0.27) 


19.00(0.30) 






15.60(0.17) 


17.21(0.19) 


30 


>18.68 


18.79(0.07) 


19.97(0.03) 


20.25(0.05) 


18.91(0.03) 






15.60(0.02) 


17.50(0.11) 


31 


>18.82 


18.97(0.14) 


20.40(0.05) 


20.65(0.25) 


19.77(0.03) 






16.58(0.12) 


18.03(0.07) 


32 


16.85(0.05) 


17.30(0.03) 


18.39(0.02) 


18.61(0.02) 


18.19(0.02) 






15.43(0.02) 


16.98(0.01) 


33 


18.86(0.38) 


19.09(0.37) 


19.75(0.10) 


20.32(0.33) 


19.09(0.07) 






16.02(0.12) 


17.68(0.17) 


34 


>18.91 


>21.19 


22.33(0.10) 


21.96(0.21) 


21.11(0.04) 






18.31(0.11) 


19.40(0.21) 


35 


>19.21 


19.43(0.07) 


20.82(0.01) 


20.98(0.16) 


20.13(0.03) 






17.09(0.19) 


18.53(0.12) 


36 


>19.08 


>21.36 


23.16(0.22) 


22.36(0.21) 


21.61(0.21) 






19.34(0.15) 


20.12(0.13) 


37 


16.91(0.14) 


17.48(0.03) 


18.98(0.01) 


18.90(0.03) 


18.49(0.00) 






15.18(0.02) 


16.96(0.04) 


38 


17.17(0.14) 


17.76(0.03) 


19.18(0.06) 


19.10(0.08) 


18.90(0.06) 






16.03(0.02) 


17.47(0.03) 


39 


17.17(0.11) 


18.15(0.05) 


19.66(0.03) 


19.78(0.05) 


19.25(0.04) 






15.87(0.03) 


17.87(0.04) 



Note. — We adopt a distance modulus of ~ 28 for NGC 4449. A three-dot ellipsis indicates no data was available. 



W^Ha = 



F, 



Hq 



„(Ha) ' 

^cont 



(1) 



where i^Ha is the total Ha flux (erg s ^ cm ^) and 
fco^t is the continuum flux density (erg s^^ cm^^ ^^^) 
at 6563 A. The F658N fllter contains both Ha and 
[NHJA6584 emission so the contribution from the [NH] 
line must be subtracted. i<Ha is given by 



f ,fF658N — 
/f660N 



/•(Ha) 
Jcont 



AA 



F658N 



(2) 



/•(NII)\ A ^ 

/cont ) ^-\F660N 



where the flrst term in square brackets is the total Ha + 
[NHJA6584 flux and the second term in square brackets 
is the total [NHJA6584 flux. Differences in throughput 
of the fllters at the wavelengths of Ha and [NH] are ac- 
counted for in Equation 2. The measured flux densities 
and widths of the filters are given by /f658N,/f660N and 
^-^F658N,AAFggoN, respectively. The flux densities of the 
continuum values, /cont and /cont ; are found by inter- 
polating between the F550M and F814W data. 



We compare the calculated Ha equivalent widths, 
VFhq, to the models in order to estimate ages for the ra- 
dio sources in NGC 4449. The thermal radio sources in 
NGC 4449 have Ha equivalent widths between ~ 340 and 
^ 2460 A, implying they have ages < 5 Myr. However, 
it should be noted that the starburst99 models reflect 
a fully sampled stellar IMF. For clusters with masses 
< 10® Mq, it becomes increasingly unlikely that a given 
cluster will contain extremely massive stars. As a re- 
sult, lower mass clusters will tend to produce less ionizing 
flux for the same age than would be predicted by simply 
scaling the starburst99 predictions. In addition, the 
possibility of contamination from an underlying or sur- 
rounding stellar population to the continuum measure- 
ments as a result of using radio-selected apertures (see 
§3.1) could artificially decrease the measured Ha equiv- 
alent widths. Therefore, the radio-detected clusters may 
be even younger than derived here (Table 5). 



4.2.2. Ionizing Fluxes 

The ionizing flux of a starburst region, Qlyci can be es- 
timated from its radio luminosity density following Con- 
don (1992): 



10 



Reines, Johnson & Goss 



Source 



TABLE 5 
Properties of the Thermal Radio Sources 



Log Who 

(A) 



Age 

(Myr) 



QLyc 
(10« s-1) 



Mass 
(103 Mq) 



Av 
(mag) 



2 


>3.10 


3 


3.00(0.28) 


4 


2.68(0.07) 


5 


2.73(0.16) 


7 


3.12(0.22) 


11 


2.53(0.50) 


18 


3.26(0.07) 


22 


2.79(0.03) 


23 


3.21(0.08) 


27 


3.13(0.14) 


29 


3.20(0.58) 


30 


3.39(0.09) 


33 


3.24(0.48) 


Note. 


— Ha equ 



<3.7 


17(3) 


(1.2)4.0(0.8) 


40(4) 


(0.2)4.9(0.1) 


39(5) 


(0.4)4.8(0.3) 


19(4) 


(0.6)3.3(1.0) 


18(3) 


(2.1)5.2(2.2) 


33(4) 


(0.2)2.9(0.1) 


21(2) 


(0.1)4.6(0.1) 


23(2) 


(0.1)2.9(0.7) 


68(3) 


(0.8)3.6(0.5) 


42(3) 


(2.9)2.9(2.1) 


23(4) 


(1.0)2.6(0.2) 


27(2) 


(2.9)2.9(1.8) 


18(3) 



<7 


4.3(0.3) 


21(2) 


0.3(0.1) 


40(5) 


0.2(0.1) 


18(3) 


1.5(0.3) 


6(1) 


1.0(0.2) 


44(5) 


1.5(0.4) 


6(1) 


0.1(0.1) 


19(2) 


0.7(0.1) 


19(1) 


1.2(0.1) 


16(1) 


0.8(0.1) 


6(1) 


1.0(0.3) 


6(1) 


1.2(0.1) 


5(1) 


1.2(0.2) 



Ha equivalent widths, ages, ionizing fluxes, masses and ex- 
tinctions of the thermal radio sources in NGC 4449. These values were 
derived from the nebular emission from the HH regions (§4.2). 



!*)s"^'°"(lli< 



GHz 



(3) 



L 



t^, thermal 



1027 erg s-i Hz-i 

Equation 3 is applicable to thermal radio emission but 
provides only a lower limit since some of the ionizing 
photons may be absorbed by dust and the radio emis- 
sion may be optically thick. Studies of Lyman contin- 
uum extinction in Hll regions find that the fraction of 
ionizing photons absorbed by dust can be quite large 
(Inoue 2001; Dopita et al. 2003), leading to derived ion- 
izing fluxes (and therefore masses) that can be artificially 
low by up to a factor of 10. Nonetheless, we adopt an 
electron temperature of Tg = 10^ K and use the 3.6 cm 
radio luminosities to calculate QLyc using Equation 3 for 
the radio sources in NGC 4449. The 1.3 cm radio lumi- 
nosities are somewhat less likely to be contaminated by 
non-thermal emission or be optically thick, but the 1.3 
cm image does not have the requisite sensitivity. The 
derived QLyc of the thermal radio sources range from 
~ 20-70 xlO''^ s-i, or equivalently (Vacca et al. 1996), 
~ 20-70 07.5 V stars. 

4.2.3. Masses 

The ionizing flux and age of a thermal radio source 
can be used estimate its total stellar mass. The models 
predict the number of ionizing photons at a given age for 
an input total stellar mass. Using the age derived from 
the Ha equivalent width, we determine the expected ion- 
izing flux at that age, given by the model. We assume 
that mass scales with ionizing flux and estimate the to- 
tal stellar mass by multiplying the input model mass by 
the ratio of the observed to expected ionizing flux. The 
masses of the thermal radio sources derived using this 
method are in the range ~ 0.5-5 xlO^ M©. 

4.2.4. Extinctions 

Global extinction of the ionized nebular gas is esti- 
mated by comparing the measured Ha flux to the pre- 
dicted Ha flux (with no extinction) derived from the 



thermal radio flux density, 5'i/^thormai- The radio emission 
does not suffer from extinction per se, although again we 
note that the UV photons can suffer significant absorp- 
tion before ionizing the gas, which is not accounted for 
in Equation 3. We follow Condon (1992) to obtain the 
predicted Ha flux: 



^ FHa,prcdictcd \ _ Q.g X 10" 

yCrg s ^ cm ^ J 



12 



Te 



-0.59 



X I 

VGHz 
The extinction of Ha is then given by 



Vio^Ky 

/ '-'l', thermal 



(4) 



mjy 



Ana = -2.5 X log 



/ Fna. 



observed 



\F}ia 



predicted 



(5) 



A 30 Doradus extinction curve (Fitzpatrick 1985; Fitz- 
patrick & Massa 1990; Misselt et al. 1999) is used to 
convert Aho to Ay, the extinction at 5500 A. The ther- 
mal radio sources have measured nebular extinctions of 
Ay ^1.5 with the exception of Source 2, which is highly 
obscured and has Ay « 4.3. Measured extinctions. Ha 
equivalent widths, ages, ionizing fluxes and masses of the 
thermal radio sources in NGC 4449 are listed in Table 5. 

4.3. Results from SED Fitting: Evidence for I- and 
H-band Excesses 

In addition to using the nebular emission, we also esti- 
mate the physical properties of the radio selected sources 
in NGC 4449 using the stellar emission by comparing a 
grid of model spectral energy distributions to the broad- 
band HST flux densities. The grid consists of models 
with ages up to 10 Myr in increments of At = 0.1 Myr 
and extinctions in increments of A Ay ~ 0.1. We typi- 
cally run the grid from Ay = to 5 magnitudes which 
yields a total of 5,151 model SEDs to compare to each 
observed SED. Galactic foreground extinction is also 
accounted for (NGC 4449 has a Galactic reddening of 
E{B -V)^ 0.019 (Schlegel et al. 1998)). 

Each model spectral energy distribution, with a given 
age and extinction, is convolved with the HST total 
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TABLE 6 
Results of SED Fitting 



Extinction Curve 


Z=0.004 


Z=0.008 




Fitting the 


UV, U, B &i 


t^Data 


30 Doradus 
Galactic 
Calzctti 


0.032(0.028) 
0.040(0.034) 
0.043(0.036) 


0.032(0.027) 
0.040(0.035) 
0.042(0.036) 



Fitting the UV, U, B, V & / Data 



30 Doradus 
Galactic 
Calzetti 



0.061(0.041) 
0.055(0.038) 
0.058(0.038) 



0.061(0.038) 
0.054(0.037) 
0.056(0.037) 



Note. — The goodness-of-fit parameter, a, and 
its standard deviation, averaged over all possible 
sources, is listed for 12 runs of SED fitting. Ex- 
cluding the /-band flux and using a 30 Doradus 
extinction curve fits the data significantly better 
than the other runs. A full discussion is given in 
the text. 



throughput curves corresponding to the broad-band fil- 
ters listed in Table 2 before it is compared to the obser- 
vations. Once we have obtained the model flux densities 
for each broad-band filter, we compute a goodness-of-fit 
statistic, a, given by the standard deviation of the differ- 
ence between the log of the observed flux densities and 
the log of the model flux densities, weighted by the mea- 
surement errors of the observed flux densities. In the case 
where we only have an upper limit for an observed flux 
density, we use this limit to exclude models but do not 
include this point in the fit. The model that minimizes 
a is chosen as the best fit. A mass estimate is obtained 
by scaling the input model mass by the weighted mean 
offset between the observed and best-fit model flux den- 
sities (again done in the log). 

We have performed SED fitting on the radio selected 
sources in NGC 4449 with detections in the HST data. 
The most heavily obscured sources are not fit since 
we only have upper limits for all of the flux densi- 
ties. The strongest Balmer lines. Ha and H/3 are not 
contained in any of the broad-band filters used here, 
although the F435W filter does include H7 and US. 
Originally, we included up to five filters in the fitting 
process: F170W (--UV), F336W (r-U), F435W (^B), 
F550M (^ V), and F814W (-^I). However, we found that 
even the best-fitting models did not visibly fit the data 
in many cases because of an apparent excess in the I- 
band. To check there were no calibration problems in 
the F814W data, we obtained WFPC2 F814W images of 
NGC 4449 from the archive. The results were consistent 
with our findings in the ACS data. 

To investigate further, we ran our SED fitting routine 
twelve times on each possible source: including and ex- 
cluding the F814W data point, using models with metal- 
licities of Z=0.004 and Z=0.008, and applying three dif- 
ferent extinction curves: the Galactic extinction curve 
(Cardelli et al. 1989), the starburst obscuration curve 
(Calzetti et al. 2000), and a 30 Doradus extinction 
curve adopted from Misselt et al. (1999, Table 3) and 
Fitzpatrick (1985, Table 6), using the parameterization 
given by Fitzpatrick & Massa (1990). For each of the 
twelve runs, we have computed the mean a (goodness- 
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Fig. 6. — Color-color diagram of the thermal radio sources in 
NGC 4449. This plot of B - V (F435W-F550M) versus V - I 
(F550M-F814W) illustrates that the reddest clusters in V - / are 
the bluest in B — V, contrary to the expected trend if extinction 
were responsible for the red V — I values (i.e. the /-band excess). 
A model evolutionary track and a reddening vector (using the 30 
Doradus extinction curve) are also shown for comparison. 



of-fit statistic) of all the sources and its standard de- 
viation (Table 6). These values measure how well a 
family of models (metallicity, extinction curve, includ- 
ing/excluding /-band data) fit the data set as a whole. 
Table 6 shows that excluding the /-band data and using 
the 30 Doradus extinction curve fits the data significantly 
(~ 20-50%) better than the other families of models, but 
that there is no significant difference between the fits us- 
ing either metallicity. 

It is also evident in Table 6 that excluding the /-band 
data significantly improves the fits, regardless of which 
extinction curve is used. A better fit is expected if using 
fewer data points (i.e., excluding the /-band); however, 
we ran Monte Carlo simulations which showed that in- 
cluding one more data point in the fits should only in- 
crease CT by ~ 2-5%, not 33-91% as seen in Table 6. 
Data was simulated by adding noise (representative of 
measurement errors) to model SEDs and then convolv- 
ing the SEDs with the appropriate throughput curves. 
The simulated data were compared to the original SEDs 
and the goodness-of-fit statistic, a, was calculated with 
and without the /-band data. 

It should be emphasized that extinction of the stellar 
continuum cannot account for the observed SEDs when 
including the /-band data. Figure 6 shows a plot oi B — V 
(F435W-F550M) versus V - I (F550M-F814W) illus- 
trating that the reddest clusters in V — I are the bluest 
in B — V , contrary to the expected trend if extinction 
were responsible for the red V — I values. The model 
evolutionary track and a reddening vector (using the 30 
Doradus extinction curve) are also shown for comparison. 
This anomalous color-color diagram, along with our re- 
sults from SED fitting, leaves little doubt that the /-band 
excess is real. 

On average, the /-band flux densities are ~ 60% higher 
(0.5 magnitudes brighter) than the predicted values from 
the best-fitting model SEDs. In the most extreme cases, 
the measured /-band flux densities exceed the model val- 
ues by more than a factor of 2.5. An even greater excess 
is present in the F160W (H) filter (also not included in 
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Fig. 7. — UV to near-IR observed flux densities (erg s~^ cm~^ A~^) and best-fitting model SEDs for the thermal radio sources. The 
blue data points {UV, U, B, V) are included in the fit but the red points (7, H) are not. Grids of a (goodness-of-fit parameter) are also 
shown for models of various ages and extinctions. The white space in the grids indicate models that were rejected based on the ultraviolet 
detection limits. A green cross indicates the best-fitting model and, for comparison, an orange cross indicates the values obtained from the 
nebular emission as described in §4.2. Source 2 is not included since it was not detected in the broad-band optical filters. 



the fits) with the data brighter than the best-fitting SED 
by an average of 1.6 magnitudes, or ~ 4x the flux. We 
will discuss possible origins of the /- and _H-band excesses 
in §5.1. 

Figure 7 shows the observed flux densities and best- 
fitting model SEDs for the thermal sources that have 



detections in the HST data. The models have a metallic- 
ity of Z=0.004, use the 30 Doradus extinction curve, and 
exclude the /- and iJ-band data in the fits. The grids 
of (J are also shown to give a sense of the uncertainty 
in the estimates of ages and extinctions. A green cross 
indicates the best-fitting model and, for comparison, an 
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orange cross indicates the values obtained from the nebu- 
lar emission as described in the previous section. Overall, 
there is good agreement between the physical properties 
derived from the nebular emission and SED fitting. 

4.4. Integrating Results from the Nebular and Stellar 
Emission 

By combining results from the nebular and stellar emis- 
sion, a self-consistent picture relating the observed prop- 
erties of young emerging massive star clusters to model 
predictions is beginning to develop. It is noted that cau- 
tion must be applied when interpreting the physical prop- 
erties of observed clusters derived from population syn- 
thesis models, since stochastic fluctuations can bias the 
results. However, we do not believe this is significantly 
affecting the results presented in this paper since the clus- 
ter masses are ^10^ M© and we are not solely relying on 
broad-band colors (Langon & Mouhcine 2000; Cervino & 
Luridiana 2004). In this section, we combine the results 
from the nebular (§4.2) and stellar (§4.3) emission and 
explore various trends in the data. 

Figure 8 shows the ratio of radio flux densities (3.6 cm) 
to the optical ( V'-band) flux densities plotted against Ha 
equivalent width. The ionizing flux, and therefore Ha 
equivalent width, are expected to decrease as a cluster 
ages. However, the ionizing flux is also a function of a 
cluster's mass whereas Ha equivalent width is roughly 
mass independent. Normalizing the radio flux densities 
(which traces the ionizing flux) by the optical flux den- 
sities removes the mass dependence. It is clear that the 
mass normalized radio flux densities and Ha equivalent 
widths do indeed decrease with cluster age. A model evo- 
lutionary track from 1-6 Myr is also shown in Figure 8. 
The model flux density at 3.6 cm was obtained from 
inserting the predicted ionizing flux from starburst99 
into Equation 3. The vertical offset between the data and 
the model evolutionary track can be attributed to extinc- 
tion in the F-band (top of Fig. 8) . The data match the 
model track extremely well after correcting for the ex- 
tinctions derived from the nebular emission (bottom of 
Fig. 8). This result suggests that the extinction to the 
stars and the gas is approximately the same, which is sig- 
nificant since the stellar light and nebular gas may have 
different geometries (i.e. the nebular emission is typically 
more extended than the stellar emission). 

Figure 9 also indicates that the extinction to the stars 
and the gas is roughly the same in most cases. The ma- 
jority of the thermal radio sources have consistent extinc- 
tion estimates for the stars and the gas, although approx- 
imately one-third of the sources have nebular extinctions 
that are ~ 1 magnitude higher than their derived stellar 
extinctions. Based on the results from Figure 8, these 
differences may simply be due to uncertainties in the ex- 
tinction estimates from SED fitting, which are largely 
dependent on the ultraviolet (and lowest quality) data. 
The SEDs also suffer a degeneracy between age and ex- 
tinction. 

However, it is worth mentioning here that discrepan- 
cies in extinction estimates of extragalactic Hll regions 
have been recognized since the 1970's. In particular, 
early Westerbork data for Hll regions in M33, MlOl and 
M51 (Israel & van der Kruit 1974; Israel et al. 1975; van 
der Kruit 1977) showed that comparing radio continuum 
emission to Ha implied about 1.2 mag excess extinction 
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Fig. 8. — Top: Ratio of the radio fiux densities (3.6 cm) to 
the optical ( y-band) fiux densities plotted against Ha equivalent 
width for the thermal sources in NGC 4449. Bottom: Same as 
above except the F-band flux densities have been corrected for ex- 
tinction using the estimates derived from the nebular gas. These 
data match the model evolutionary track extremely well, illustrat- 
ing that ionizing flux, and therefore H« equivalent width, do indeed 
decrease as a cluster ages (as predicted) . The model flux density at 
3.6 cm was obtained from inserting the predicted ionizing flux from 
STARBURST99 into Equation 3. These plots also suggest that the 
extinction is approximately the same to the stars and the gas. 

as compared to the Balmer decrement procedure. These 
results combined with other determinations were sum- 
marized by Israel & Kennicutt (1980), who discuss two 
possible origins for the "missing dust": 1) the dust may 
be internal and mixed with the ionized gas, 2) over large 
sizes of many hundreds of parsecs, a non-uniform dis- 
tribution of gas and dust will cause the Balmer decre- 
ment to be mainly weighted towards regions with lower 
extinction, while the radio derived absorption will be de- 
termined mainly by the more absorbed regions. It is 
unclear whether the discrepent extinctions derived for 4 
of the thermal radio sources in this work are a result of 
these types of effects, or simply due to uncertainties in 
the SED derived extinctions. 
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Fig. 9. — Comparison of the extinction estimates of the ther- 
mal radio sources in NGC 4449 derived from the nebular emission 
and SED fitting. A line showing the one-to-one correlation is also 
shown, illustrating that many sources have similar extinction es- 
timates, yet some sources have a derived nebular extinction that 
is ~ 1 magnitude higher than the stellar extinction. This discrep- 
ancy may simply be due to uncertainties in the estimates from SED 
fitting. 
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Fig. 10. — Comparison of the age estimates from H« equivalent 
widths and SED fitting (and the line showing a one-to-one correla- 
tion). Ha equivalent width is a more sensitive age indicator, since 
ionizing flux is a strong function of time once the most massive 
stars begin to die. The SEDs, on the other hand, do not change 
very much at these young ages which can explain the pile-up at 
~ 3 Myr. 



The age estimates from SED fitting also have rela- 
tively high uncertainties (see Fig. 7), and tend to be 
concentrated around ~ 3 Myr. These results are proba- 
bly due to an insensitivity of the model SEDs, which do 
not change significantly at these young ages. Ha equiva- 
lent width, however, is a more sensitive age indicator at 
young ages since ionizing flux is a strong function of time 
once the most massive stars begin to die (although these 
age estimates may be systematically high for reasons dis- 
cussed in §4.2.1). A comparison of the age estimates from 
SED fitting and Ha equivalent widths is shown in Fig- 
ure 10. 

We find two interesting correlations with Ha equivalent 
width, and hence cluster age. First, the /-band excess 
(§4.3) is largest in the youngest clusters and decreases 
with age. We define the /-band excess in magnitudes 

as A/ = -2.5 X log(/<$bLvcd//midei)> where /.^bLvcd is 
the observed flux density in the F814W filter and j\aa^^i 
is the flux density of the best fitting model SED. Fig- 
ure 11 shows a plot of A/ versus age for the thermal 
radio sources and the best linear fit is given by 
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Setting A/ equal to zero implies that whatever is caus- 
ing the /-band excess in these young clusters is unlikely 
to affect clusters older than ~ 5 Myr. We will discuss 
possible origins of the /-band excess in §5.1. 

We also observe a puzzling trend that seems to indi- 
cate cluster mass is positively correlated with age for 
ages < 5 Myr. Figure 12 shows a plot of the extinction 
corrected My versus Ha equivalent width. The model 
evolutionary tracks for various cluster masses are also 
shown. A strong anti-correlation between My and Ha 
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Fig. 11. — Difference between the observed and model /-band 
magnitudes versus age, demonstrating that younger sources have 
a larger /-band excess. The best linear fit to the thermal sources 
is shown and given by Equation 6. Based on these data, we do not 
expect the cause of the observed /-band excess to affect clusters 
older than ^ 5 Myr. 



equivalent width is evident, implying more massive clus- 
ters tend to be older than less massive clusters. The same 
overall trend is observed without correcting My for ex- 
tinction, while a positive correlation is found between 
My and Ha total flux (as expected). An observational 
bias could possibly explain why optically faint clusters 
with low Ha equivalent widths are not detected, since 
these clusters may be less likely to show up in the radio 
and make it into our sample. However, it is curious that 
we do not detect clusters that are relatively bright in V 
with high Ha equivalent widths. It is not clear what is 
causing this trend in the data, but we explore various 
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Fig. 12. — Top: Extinction corrected absolute V magnitude 
versus Ha equivalent width of the thermal radio sources. The 
extinction estimates come from comparing the radio and Ha emis- 
sion from the ionized gas. Model evolutionary tracks for a range 
of masses are also shown. This plot suggests there is a correlation 
between cluster mass and age, in excess of the model predictions. 
See §5.2 for a discussion. Middle: The same overall trend is ob- 
served without correcting M^^ for extinction. Bottom: yiy and Ha 
luminosity arc positively correlated, as expected. 

possibilities in §5.2. 
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Fig. 13. — Comparison of the stellar masses of the thermal 
radio sources in NGC 4449 derived from nebular emission and SED 
fitting. The line showing a one-to-one correlation is also shown. 
A slight trend is apparent that suggests mass estimates from the 
nebular emission tend to be higher than those from SED fitting for 
the lower mass clusters, while the mass estimates from SED fitting 
tend to be higher than those derived from the nebular emission for 
the higher mass clusters. 

Overall, the mass estimates derived from the nebular 
emission and SED fitting are consistent within a factor 
of ~ 3. However, a slight trend is evident in Figure 13 
suggesting that mass estimates from the nebular emission 
tend to be higher than those from SED fitting for the 
lower mass clusters, while the mass estimates from SED 
fitting tend to be higher than those derived from the 
nebular emission for the higher mass clusters. We will 
return to this in §5.2 

5. DISCUSSION 
5.1. Possible Origins of the I-band Excess 

We have demonstrated that the /-band fluxes of the ra- 
dio selected clusters in NGC 4449 exhibit an excess com- 
pared to model spectral energy distributions (Fig. 7). In 
the B — V versus V — I color-color diagram, the reddest 
clusters in V — I are simultaneously the bluest in B — V 
(Fig. 6). The /-band excess is also positively correlated 
with Ha equivalent width, implying the excess decreases 
with age (Fig. 11). Based on these observations, we do 
not expect the origin of this /-band excess to affect clus- 
ters older than ~ 5 Myr (Eqn. 6). 

Possible origins of the /-band excess are discussed be- 
low, but first we note that other studies have also re- 
vealed the presence of an /-band excess in a small number 
of young massive clusters. Whitmore & Schweizer (1995) 
discovered about a dozen extremely red objects in the SE 
quadrant of the merging Antennae galaxies. This region 
also contains the strongest radio continuum sources in 
the system as well as strong mid-IR emission and a few 
of these sources are probably associated with the red ob- 
jects. Whitmore Sz Schweizer (1995) interpreted these 
red objects as extremely young, embedded clusters that 
have not yet cleared out their natal dust cocoons. 

In their study of the circumnuclear starburst galaxies 
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NGC 5248 and NGC 1512, Maoz et al. (2001) found a 
number of young massive clusters with excess emission in 
the _ff-band, and one of these sources had a detectable I- 
band excess. Maoz et al. (2001) were able to fit the SED 
of this source by adding a 2000 K blackbody to their 
model and determined that if the emission were due to 
dust heated by stars, it must be circumstellar. 

Possible origins of the observed /-band excess in the 
young radio-detected clusters of NGC 4449 are consid- 
ered here. 

i) Red Supergiants 

One possible origin for an /-band excess is the presence 
of red supergiants. These evolved stars are known to 
dominate the near-infrared light of massive clusters older 
than ~ 7 Myr (e.g., McCrady ct al. 2003), but we cannot 
definitively rule out their presence in the younger clusters 
of NGC 4449 because of stochastic variations. 

In addition, since we use the 3.6 cm radio image as 
our reference when doing photometry, the sizes of our 
apertures are larger than the linear resolution of the HST 
images and individual compact star clusters. Some of 
our apertures contain multiple discrete sources, of which 
some are very red and could be individual supergiants. 

On the other hand, there are radio sources with single 
dominant optical counterparts. For example. Source 30 
is a compact optical cluster with a linear size ^ 3 pc 
(FWHM), and it exhibits the largest /-band excess (and 
youngest age). It would be quite unexpected for this 
young (< 3 Myr) cluster to contain red supergiants. In 
addition, the /-band excess decreases with cluster age 
(Fig. 11) while the presence of supergiants should in- 
crease with time for clusters less than 10 Myr old. We 
therefore consider it highly unlikely that red supergiants 
are the primary cause of the /-band excess. 

ii) Emission Lines 

Emission lines are another possible source of the ob- 
served /-band excess. The Paschen series (9 and higher) 
and the [SIII] A9069, A9532 emission lines could con- 
tribute to the fiux in the F814W filter and these lines 
have all been observed in the spectra of Hll galaxies and 
starbursts (e.g. Kehrig et al. 2006; Marquart et al. 2007). 

To see if emission lines could plausibly be the dom- 
inant source of the /-band excess, we calculate the ex- 
tinction corrected flux ratio, relative to Ha, needed to 
account for the excess flux in the /-band and find it to 
be ~ 0.55 on average with a range between ~ 0.1-1. 
In other words, the requisite contribution from emission 
lines in the F814W filter would need to be - 10-100% of 
the Ha flux for a given source. The higher level Paschen 
lines are not expected to contribute nearly this much flux 
(e.g. Pa9/Ha < 1%). The ratio of the [SIII] lines to Ha 
has been observed to be as high as 20-30% (Kehrig et al. 
2006) and may play some role here, but the integrated 
system throughput for the F814W filter drops by 50-75% 
from its maximum at the wavelengths of these lines. For 
all of these reasons, we believe that emission lines are 
probably not the dominant source of the /-band excess. 

iii) Thermal Emission from Hot Dust 

It is possible that thermal emission from hot dust 
near the sublimation temperature (^ 1700 K) could con- 
tribute to an excess in the /-band. An excess at ^ 2 /im 



is observed in many star forming regions and usually at- 
tributed to hot dust. However, the youngest radio de- 
tected clusters in NGC 4449 have such large /-band ex- 
cesses that it would be truly remarkable if thermal emis- 
sion from heated dust were the dominant cause since 
either there must enough emitting dust near the sub- 
limation temperature to dominate the stellar light at 
~ 0.8 /xm or the dust must survive even higher tem- 
peratures. 

iv) Optically Veiled Photospheric Continuum from Em- 
bedded Stars 

A sub-population of deeply embedded stars that have 
not yet fully emerged from their protostellar clouds could 
also potentially cause an /-band excess. These stars 
would be invisible at optical wavelengths because of high 
extinction, but their photospheric emission could con- 
tribute to the integrated cluster light in the near-IR and 
even in the /-band. However, the observed /-band ex- 
cess is so large in some sources with no obvious excess 
in the F-band, that the embedded population would be 
required to simultaneously have moderate to large ex- 
tinctions (^v ~ 5) and uncomfortably high fractions 
(^ 90%) of the total cluster mass. 

v) "Extended Red Emission" 

Extended Red Emission (ERE) is a very promising ex- 
planation for the /-band excess (see Witt & Vijh (2004) 
for a review) . ERE is typically observed as a broad emis- 
sion feature (^ 1000 A) which peaks in the red to near-IR 
spectral range. ERE has been discovered in number of 
astrophysical environments where both interstellar dust 
and ultraviolet photons are present, including the Orion 
Nebula (Perrin & Sivan 1992) and the 30 Doradus Neb- 
ula (Darbon et al. 1998). 

ERE is caused by photolumincscense, a 3-step process 
in which 1) an interstellar particle absorbs an ultraviolet 
photon and makes an electronic transition to an excited 
state, 2) a series of rotational/vibrational transitions re- 
laxes the system to an intermediate state, and 3) an op- 
tical photon is emitted through an electronic transition 
back to the ground state. The carrier of ERE is unknown. 

Studies of ERE have revealed that, in any given envi- 
ronment, the maximum intensity is well correlated with 
the local density of the radiation field (Gordon et al. 
1998), and in particular with the density of the ultra- 
violet radiation field (Smith & Witt 2002). The ERE 
peak wavelength is also closely correlated with the den- 
sity of far-UV photons, suggesting that the radiation 
field is not only driving the ERE, but also modifying 
the carrier through size-dependent photo-fragmentation, 
preferentially affecting the smallest carrier particles first 
(Smith & Witt 2002). Furthermore, nebulae in which 
ERE is not detected (Witt & Boroson 1990) tend to 
be diffuse (unclumped) and have had exposure to UV 
radiation for an extended amount of time. In environ- 
ments with clumpy dust, ERE is strongest on the faces 
of high-density clumps and undetectable in low-density 
interclump regions. 

These findings are consistent with ERE causing the 
observed /-band excess in the young massive star clus- 
ters in NGC 4449 (recall that the excess decreases with 
cluster age and is expected to disappear after ~ 5 Myr). 
The youngest clusters will have the most intense UV ra- 
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diation fields and they are also the most likely to have 
a clumpy dust structure, both of which would increase 
the intensity of the ERE and hence the /-band excess. 
Over time, the clumpy dust structure may be dispersed 
and/or the radiation field may modify/destroy the ERE 
carrier, both causing the /-band excess to decrease and 
then disappear all together. 

ERE offers a favorable hypothesis for the /-band ex- 
cess, although spectroscopic follow-up is required for con- 
firmation. It is also important to note that if ERE is the 
cause of the /-band excess, the //-band excess likely has 
a different origin. Contributions from hot dust, embed- 
ded stars and red supergiants are all possible. We are 
currently undertaking follow-up spectroscopic and pho- 
tometric observations of NGC 4449 to further investigate 
the /- and //-band excesses. 

5.2. Interpretation of the Mass-Age Correlation for 
Ages < 5 Myr 

Figure 12 illustrates that while My and Ha flux are 
positively correlated (as expected). My and Ha equiv- 
alent width are anti-correlated. In other words, fainter 
clusters have relatively more Ha emission with respect 
to the stellar continuum than brighter clusters. This can 
be interpreted as a mass-age correlation, since My gen- 
erally traces cluster mass (modulo extinction) and Ha 
equivalent width traces age. The two most trivial, and 
unsatisfying, interpretations of this trend are 1) the low- 
mass clusters formed more recently than the high-mass 
clusters, and 2) smaller clusters coalesce and make larger 
clusters over time. However, we note that there is no 
observed correlation between cluster size (measured by 
aperture size) and Ha equivalent width. 

A third interpretation of the apparent mass-age cor- 
relation for clusters in this work (with ages < 5 Myr), 
although purely speculative, is the following. Is is possi- 
ble that the observed correlation is actually the imprint 
of an evolutionary sequence, where the youngest clus- 
ters merely appear less massive than their older counter- 
parts. First, relative to more evolved clusters, younger 
clusters are likely to contain a higher fraction of dust en- 
shrouded, optically veiled stars. These heavily embedded 
stars would not contribute to the integrated F-band flux, 
yet the young optically visible stars within the clusters 
would still produce large Ha equivalent widths revealing 
their young ages. Therefore, younger clusters may tend 
to appear less massive. In addition, if the younger clus- 
ters have a majority of their mass in the optically veiled 
sub-population of stars, it would not be surprising that 
the masses inferred from the nebular emission (i.e. radio 
and Ha equivalent width, which are mostly unaffected 
by extinction) are higher than those inferred from the 
stellar emission (which only samples the optically visible 
"emerged" stars) as is apparent in Figure 13. 

A second factor that could contribute to the appar- 
ent mass-age correlation is stellar winds. More massive 
clusters will have disproportionately strong stellar winds, 
which could preferentially clear away the gas. A lower 
gas fraction will result in a lower equivalent width, and 
therefore a larger inferred age. As a result, more massive 
clusters would tend to appear older by virtue of their 
dearth of ionized gas. A deficiency of ionized gas could 
also explain why the more massive clusters tend to have 
higher mass estimates from SED fitting than those de- 




FlG. 14. — VLA 3.6 cm radio image of NGC 4449. A 1.6 kpc su- 
perbubble may iiave triggered the formation of tiic 7 radio sources 
at tile circle's perimeter. See §5.3 for more discussion. 

rived from the nebular emission (Fig. 13). 

5.3. Evidence for Triggered Star Formation in 
NGC 4449 

NGC 4449 is thought to have undergone some type 
of interaction in its past (Hunter 1998; Theis & Kohle 
2001), which has likely influenced the global star forma- 
tion history of the galaxy over the last ~ 1 Gyr (e.g. An- 
nibali et al. 2007). Local processes such as stellar winds 
and supernovae (perhaps the result of a previous burst) , 
however, seem to be playing a significant role in trigger- 
ing the current bursts of star formation, resulting in the 
formation of massive clusters detected in the radio. 

Figure 14 illustrates that the formation of a number of 
radio sources appears to have been triggered by an ex- 
panding superbubble. Taken alone, the fact that seven 
sources (23, 24, 25, 28, 30, 37, & 39) he along the same 
circular arc (radius = 44"= 832 kpc) could be viewed as 
a coincidence, albeit an unlikely one. However, previous 
observations of gas and dust in NGC 4449 suggest oth- 
erwise. Hi, CO and cold dust appear to be depressed or 
almost absent in this region of the galaxy (e.g. Fig. 1 
in Bottner et al. 2003) while soft diffuse X-ray emission 
fills it in (Fig. 5 in della Ceca et al. (1997) and Fig. 10 in 
Summers et al. (2003)). In addition, kinematic studies of 
NGC 4449 are consistent with the stellar component of 
the galaxy observed face-on (Hunter et al. 2002, 2005), 
which would account for the sources lying on a circular 
arc in projection. All of this evidence supports the view 
of an expanding superbubble triggering several sites of 
star formation in NGC 4449. 

Further evidence of triggered star formation is appar- 
ent within NGC 4449. An Ha shell with radius ~ 250 pc 
is observed in the SW region of the galaxy (see Figs. 3 & 
4) surrounding numerous optically visible stars and clus- 
ters, whose stellar winds and supernovae likely produced 
the shell, while six radio sources (1, 2, 3, 4, 5, & 7) are 
found along the periphery of this shell. It is possible that 
a collision between this shell and the larger superbubble 
is responsible for compressing the gas in the southern re- 
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gion of the galaxy and triggering the formation of Sources 
18, 20, and 21. 

Another group of radio sources (31-36) is found along 
a smaller Ha loop (radius ~ 40 pc) in the NE region 
of the galaxy, again containing optically visible stellar 
clusters that could have triggered the current local burst 
of star formation. 

5.4. Radio Sources without Optical or Infrared 
Counterparts 

A close inspection of Figures 3 & 4 reveals that at least 
a few radio sources in NGC 4449 do not have optical 
counterparts (e.g. Sources 20, 24, & 28 ). Sources 24 & 
28 also do not appear to have any detectable 24 /im emis- 
sion associated with them. It is unclear whether this is 
the case for the other radio sources without optical coun- 
terparts due to confusion in the 24 (mii image. In addi- 
tion, the "pure" radio sources have only been detected 
at 3.6 cm and we cannot determine whether they are 
thermal. We do believe these sources are in the galaxy 
and not background objects because their locations are 
associated with star forming regions. With the current 
data, however, we can only speculate about the nature 
of these sources. 

If the sources are thermal, the lack of Ha could be ex- 
plained by heavy extinction from dust. For Sources 24 & 
28, the extinctions would have to be at least ^ 7 visual 
magnitudes (§4.2.4). Again, if we assume these sources 
are thermal, they are probably younger than the radio 
detected clusters with optical counterparts (< 2.5 Myr). 
They would also have relatively low ionizing fluxes (the 
equivalent of ~ 6-15 07.5 V stars) and masses. These 
young ages and low masses would be consistent with the 
observed correlation between derived masses and ages 
for the thermal radio sources with optical counterparts. 
However, if the visibly obscured sources are heavily em- 
bedded clusters, we would expect them to produce ther- 
mal emission from heated dust and be observable in the 
mid-IR which does not seem to be the case (at least for 
Sources 24 & 28). It may be more likely that these "pure" 
radio sources are supernova remnants, since non-thermal 
emission is anti-correlated with bright mid-IR emission 
(e.g. Brogan et al. 2006). 

6. SUMMARY AND CONCLUSIONS 

We have presented a multi- wavelength study of emerg- 
ing massive star clusters in the irregular starburst galaxy 
NGC 4449. We combine sensitive, high resolution 
VLA radio observations with ultraviolet, optical and in- 
frared data from the Hubble and Spitzer Space Telescope 
archives to obtain a comprehensive view of star formation 
occurring within the galaxy and to determine the phys- 
ical properties of the radio-detected star clusters. We 
summarize the results of our study below: 

1. We identify 13 thermal radio sources in NGC 4449. 
These sources have ages < 5 Myr (based on 
Ha equivalent widths), ionizing fluxes between 
~ 20-70 xlO'^^ s^i (the equivalent of - 20-70 
07.5 V stars), and stellar masses between ^ 0.5- 
5 xlO* M0. 

2. Of the 13 thermal radio sources, 12 have measured 
5/0 &aZ extinctions < 1.5 magnitudes at 5500 A and 



the extinction appears to be approximately the 
same to the stars and the gas. The most optically 
obscured thermal radio source has an Ay w 4.3. 

3. Model SEDs modified by a 30 Doradus extinc- 
tion curve fit the clusters' broad-band ultravio- 
let and optical flux densities {UV,U,B,V) signif- 
icantly better than the standard Galactic extinc- 
tion curve and the starburst obscuration curve of 
Calzetti et al. (2000). 

4. The clusters exhibit an /-band excess (with respect 
to model SEDs) that is anti-correlated with age. 
In the most extreme cases, the /-band magnitudes 
are more than 1 magnitude brighter the best-fltting 
model SEDs. We do not expect clusters with ages 
> 5 Myr to have this excess. A photoluminescent 
process, known as Extended Red Emission, pro- 
vides a favorable hypothesis for the cause of the 
observed /-band excess. 

5. The 7 nuclear clusters with NICMOS data exhibit 
a large (~ 1.6 mag) //-band excess. We speculate 
that some combination of red supergiants, thermal 
emission from hot dust, and optically veiled photo- 
spheric emission are the most likely causes. 

6. An apparent mass-age correlation (for ages < 
5 Myr) is observed for the thermal radio sources 
in NGC 4449. 

7. The mid-infrared (24 /im) morphology of 
NGC 4449 is almost identical to that of the 
radio, indicating that nearly all of the warm dust 
in the galaxy is associated with the most recent 
episodes of star formation in the galaxy. 

8. Local processes such as supernovae and stellar 
winds are likely playing an important role in trig- 
gering the current bursts of star formation within 
NGC 4449. 

Synthesizing high quality data sets across the electro- 
magnetic spectrum has provided us with a unique and de- 
tailed view of the current star formation occurring within 
NGC 4449. Our findings would not have been possible 
without the combination of radio, optical, and infrared 
observations. The results presented here are the first of 
a larger program to study massive star clusters as they 
emerge from their birth material and transition from be- 
ing visible in the radio to optical wavelength regimes. 
Ultimately, understanding the formation and early evo- 
lution of massive star clusters will shed light on a funda- 
mental mode of star formation throughout the universe 
and provide insight into the origin of globular clusters. 
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APPENDIX 
SURPHOT: A MULTI- WAVELENGTH PHOTOMETRY CODE 

We have developed a new multi-wavelength photometry code (SURPHOT) to measure flux densities within identical 
irregular apertures in multiple images. SURPHOT is written in the Interactive Data Language (IDL) and the basic 
algorithm is as follows: The user supplies a reference image and selects an aperture type, including "contour" or "free- 
form" (in addition to a circle or an ellipse). If a contour aperture is chosen, SURPHOT will detect a specified contour 
level to define a region of interest. Alternatively, the user can choose to draw a free- form aperture by clicking around a 
region of interest. Once an aperture is defined in the reference image, a background annulus is assigned by expanding 
the aperture out radially by two specified factors, matching the overall shape of the aperture. The pixel coordinates of 
vertices in the aperture and background annulus are transformed into right ascension and declination using astrometry 
keywords in the FITS header of the reference image. The RA and DEC values can then be transformed into pixel 
coordinates of any other image using its own FITS header astrometry keywords. In this way, SURPHOT uses the 
same aperture and background annulus to compute fiux densities in any number of user supplied images, provided 
they are registered to match the astrometric reference frame of the reference image. 

For a given source, 4 background levels and their errors are estimated in each image from which the user can 
choose: mode, mean, resistant mean and low mean. The mode and resistant mean are calculated using the MMM 
and RESISTANT_MEAN routines, respectively, which can be found in the IDL Astronomy User's Library ^. The 
low mean is found by using MEANCLIP (also found in the IDL Astronomy User's Library) to compute a 3a clipped 
mean of the 25% lowest valued pixels in the background annulus. In general, the mode is our preferred estimate of the 
background level, however, the other estimates can provide a consistency check or an alternative value if appropriate. 
The estimated background value (per pixel) is multiplied by the number of pixels in the aperture and then subtracted 
from the sum of the pixels in the aperture to obtain the total background subtracted source counts. In the case of the 
VLA images, the image units of Jy beam~^ are first converted to Jy pix~^, where the conversion factor is given by the 
pixel area of the 2-D Gaussian beam. The background subtracted fiux density (in Jy) is then computed as described 
above. 

Error estimates include contributions from three sources: Poisson noise from the source counts, Poisson noise from 
the background counts, and the error in measuring the background level. These error terms are calculated in units of 
electrons, added in quadrature and converted back to counts using the gain. A Icr detection limit is found in the same 
way, but only includes contributions from the Poisson noise from the background counts and the error in measuring 
the background level. 
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